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Topological effects on transition temperatures and response functions in 
three-dimensional Fermi superfluids 
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We investigate the effects of topological order on the transition temperature, T^, and response 
functions in fermionic superfluids with Rashba spin-orbit coupling and a transverse Zeeman field in 
three dimensions. Our calculations, relevant to the ultracold atomic Fermi gases, include fluctuations 
beyond mean-field theory and are compatible with /-sum rules. Reminiscent of the px + ipy super¬ 
fluid, the topological phase is stabilized when driven away from the Bose-Einstein condensation and 
towards the BCS limit. Accordingly, while experimentally accessible, is significantly suppressed 
in a topological superfluid. Above Tc, the spin and density response functions provide signatures of 
topological phases via the recombination or amplification of frequency dependent peaks. 


Introduction.— The excitement surrounding topologi¬ 
cal superfluids w derives from both their scientific as 
well as technological potential. Inspired by the canonical 
topological superfluid, a spinless Px + ipy superfluid [T], 
it has been argued m that some combination of spin- 
orbit coupling (SOC), Zeeman field, as well as superfluid 
pairing can artificially produce such a state. This was 
explored via the proximity effect in solid state sys¬ 
tems and using intrinsic pairing in ultracold atomic Fermi 
gases [HllSHTT]. There is, moreover, widespread interest in 
experimental confirmations of topological band-structure 
signatures [IMS]. 

This leads to the goals of this paper. Here we aim to 
determine (in the case of intrinsic pairing) how a tran¬ 
sition from a trivial to a topological phase is reflected 
in the superfluid transition temperature T^. Addition¬ 
ally, we show how and when a transition in the band- 
structure can be experimentally detected via studies of 
the finite temperature density-density and spin-spin cor¬ 
relation functions. The ultracold atomic Fermi gases are 
ideally suited for tuning between trivial and topological 
phases. As in the Px + ipy superfluid, from the perspec¬ 
tive of Tc, we find these intrinsically paired superfluids 
self-consistently adjust to stabilize topological phases in 
the BCS regime. This occurs despite the fact that mod¬ 
erate SOC reinforces BEC behavior through enhanced 
pairing [nMsi, even in the normal state [HIlIlllH. 

Experiments require the consideration of non-zero tem¬ 
perature T. Although studies of the ground state have 
been the focus [T2| , finite T effects have been included in 
the literature [UllSj at the mean-field level. A major 
weakness of this approach is that computing Tc in this 
manner does not reflect the topological band-structure, 
which depends importantly on the existence of a nor¬ 
mal state pairing gap. Here we remedy this inadequacy 
through the inclusion of fluctuations HU HD, and also 
establish that Tc is experimentally accessible. 

There are proposals in the literature which suggest that 
the topological phase might be observed in atomic Fermi 
gases through the compressibility k mm or via radio 
frequency (RF) based probes [H]. However, changes in 


K appear to reflect topology only in the limit of small 
SOC [131 E]. RF experiments in principle measure the 
electronic dispersion, but resolution and finite tempera¬ 
ture broadening effects are not yet HS] well controlled. 
Here we introduce an alternative probe: the frequency 
dependent density-density or spin-spin correlation func¬ 
tions |26j at temperatures T > Tc- The position or 
threshold of peaks in these responses, importantly, re¬ 
flects band-structure. In the topological phase we find 
that a peak in the density response is significantly ampli¬ 
fied due to a saddle point Van Hove singularity, often seen 
in correlated superfluids HZlHHj. In the trivial phase the 
spin response exhibits two distinct peaks, which merge 
into a single peak in the topological phase. 

Background theory.— We consider a Fermi superfluid 
described by the single particle Hamiltonian iLo(k) = 
^k+h (k).(T, where ^k = describes a free parti¬ 

cle of momentum k = (kx, ky, kz), mass m, and chemical 
potential p. The vector h (k) = hj_ (k) -|- h|[ (k) couples 
the spin-1/2 operator cr — {<Jx, Oy, Uz) to a Zeeman field 
via h|[ (k) = bzZ and an in-plane SOC field h_L (k) = 
Akj^/m for in-plane momentum kj^ = (kx,ky,Q) and 
SOC strength A. Throughout we set li = fcs = 1- 

The many-body Hamiltonian is of the Bogoliubov-de 
Gennes (BdG) form 

_^^(k)), (1) 

where A is a pairing gap and Hq (k) = (jy [Hq (— k)] ay 
is the time-reversed single-particle (hole) Hamiltonian. 
There are four branches in the BdG eigenvalue spectrum, 
vEak for a,r] = ±1 with the positive energy dispersion 

Eo.k = + + ( 2 ) 

A three-dimensional superfluid described by the above 
BdG Hamiltonian belongs to one of three distinct topo¬ 
logical phases. The topological phase diagram is spec¬ 
ified by inequalities derived from solving E-{kz,k± = 
0) = 0. No nodes appear when bz < A, corresponding 
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to a non-topological or “trivial” superfluid. If /x > 0 
and (/x^ + A^) > > A^, the topological super¬ 

fluid has four nodes (4-Weyl points) which emerge at 
kl = fx ± — A^. Finally, for arbitrary /x, the sys¬ 

tem is a topological superfluid with two nodes (2-Weyl 
points) when bl > {y? + A^) [T31 [HI E5] . For Rashba 
SOC, the dispersion around the nodes is linear in mo¬ 
mentum, and is described by a Weyl Hamiltonian with 
topologically protected nodes. 

To compute the transition temperature T^, we build 
on the well established mean-field theory miiiiiia, and 
incorporate fluctuation effects in a consistent fashion [Ml 
I29j . We write the mean-field gap equation [nKiaEs] as 

p-l j-Qj _1 ^ ^ — ivf (Egw) + / (^g'k)) \ 

2 , , V “1“ Ca'k / 

k r}OLOc' 

( 1 ^; 1 ^) 9 — (^) 

where f{x) is a Fermi distribution and g < 0 is an at¬ 
tractive interaction. Where relevant, we regularize inte¬ 
grals by introducing a scattering length defined through 
= m/47ra — [SO]- The coherence factor 

Vrjaa' (k, k) (and its generalization, Vjjaa' (k,k —q)), is 
presented in the Supplemental Material m- Their spe¬ 
cific form is irrelevant for the present discussion. 

One has to distinguish Tc from the lowest tempera¬ 
ture, denoted T*, at which the mean-field gap equation 
satisfies A(T*) = 0. Such an analysis requires a natural 
extension pn [M] of Eq. (§ to finite Q = {iuj, q) (where 
iu! is a Matsubara frequency): 

p-l (Q) ^ ~ (r;/(E„0 + /(g,,i,,q)) \ 

2 “ V + ^a'k-q) — / 

k ijaa' ^ • = -1/ / 

X v,jaa> (k,k - q) + g~^ . (4) 

From the structure of Eq. ^ it is apparent that F (0) 
depends on both the full energy spectrum Eak as well 
as the bare energy ^^k- Thus, one might expect (as im¬ 
plemented in Eq. 0) , that the fluctuation corrections 
should depend on an asymmetric combination of bare 
and dressed Green’s functions |5^ . 

The quantity F {Q) has been of interest [T71 |M] for 
computing the binding energy and mass of the pairs asso¬ 
ciated with the isolated two-body physics. We emphasize 
that to describe pairs which are intrinsic to the many- 
body system, one should not set A and fx to zero, as is 
usually done dZlIM]; the many-body state is not simply 
a gas constructed from entities of the two-body problem. 
In 3D there are metastable or resonant pairs for all pa¬ 
rameters below T*, whereas the two-body bound states 
only exist for positive scattering length. Furthermore, for 
sufficiently large b^, the effective mass in the two-body 
problem diverges near unitarity [23], and is not defined 
at negative scattering length. 

To characterize the mass of the many-body system, 
consider the vertex in E q. 0 expanded at small mo¬ 
menta, where (using Eq. 0) F (Q) ~ Oq ^(xw — Wq)”^, 


with oq = (^icjF ((5 )|q^p. The pair dispersion is Wq = 
Oq ^ (F (0,q) - F (0,0)) « q\/2Mj_ +qy2M\\, where M±_ 
(M||) is the effective pair mass for the component of 
momentum parallel (perpendicular) to the SOC vector. 
While it is sometimes possible to calculate the effective 
pair masses Mj_, My analytically, in general this is not 
necessary. Rather, it suffices to calculate numerically the 
second-order derivative at small q [33] . 

Notice that the vertex function F {Q) reflects, up to a 
constant, the non-interacting Green’s function of a ther¬ 
mal Bose gas with pair dispersion Wq, below the con¬ 
densation temperature (F”^ (0) = 0). We interpret this 
Q 7 ^ 0 contribution by considering the following quantity, 
proportional to a Bose occupation number ub'- 

^ F (Q) = A2 = ay ns, at T = T^. (5) 

Q#o 

This characterizes the excitation gap in the limit in 
which all pairs are non-condensed, but in which Tc is ap¬ 
proached from below. The Thouless condition (describ¬ 
ing the instability of the normal state) or divergence of 
F(0) requires that the above expression for A^ lies on the 
mean-field curve determined from Eq. 0 . The condition 
for Tc is then simply obtained [2Tl[Mll2^ by equating the 
constraint on A{Tc) via Eq. 0 with that obtained from 
the mean-field gap equation in Eq. 0. 

Importantly, all these arguments can be generalized so 
that the computation of Tc as a fluctuation correction 
to any BGS-BEC mean-field theory (including LOFF- 
like phases) is now accessible, based on writing F((5) as 
a natural extension of the appropriate mean-field vertex 
function (e.g., Eq. 0). This approach is distinguished 
from other BGS-BEC theories [MMS] by the fact that the 
transition temperature is dependent on a normal state ex¬ 
citation gap; in this way Tc will reflect changes in band- 
structure associated with the transition from trivial to 
topological phases. It is also distinguished by the fact 
that the present approach avoids the unphysical first or¬ 
der transition found in all other BGS-BEC theories m- 
We will assume throughout that, above Tc, the mean- 
field gap represents a reasonable approximation [38] for 
the normal state A. 

Phase diagram.— To understand the effects of SOC 
and the Zeeman field on condensation and pairing, 
we numerically compute T* and Tc, varying Xjkpa, 
bz, and A. Where relevant, we measure quantities 
in terms of the Fermi momentum {kp), energy {Ep), 
or temperature {Tp). It is convenient to define a 
shifted chemical potential, d/x = /x — /xq, where go = 
-ma.x{Eso{l + bl/E^c)/‘^ybz} and E^o = j m is the 
SOC energy. In this way, a necessary (but not sufficient) 
condition for a topological phase is that Sg > 0. 

Figure[^plots Tc and (in some cases) the pairing onset 
temperature T* as a function of either A, b^, or 1/kpa. 
Dotted lines indicate where 6g = 0. Where relevant. 
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(a) Phase diagram {b^Ep=2.5, l/kpa=0) 
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(c) Phase diagram (Ay^:^l, l/kpa=-l) 



(b) Phase diagram (b^Ep=0.5, k/kj—O.S) 
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(d) Phase diagram (kJkj:=l, l/kpa=2) 
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FIG. 1. Phase diagrams for the superfiuid temperature Tc 
and (where shown) the mean-field transition temperature T*. 
The T = 0 topological phases are indicated by shaded regions 
in light (dark) blue color with 2 (4)-Weyl points. The top two 
panels show the dependence on either SOC strength A (left) or 
scattering length l/kpa (right), with other parameters fixed. 
The lower panels show Tc vs. 6z with A/A:_f = 1 for both 
panels, Ijkpa = —1 on the left and l/kpa = 2 on the right. 
Dotted lines indicate 5^ = 0. Once a topological phase is 
entered the system becomes more BCS-like. 


these plots are consistent with earlier work El- A close 
analogy between varying A and varying l/kpa is seen in 
Fig. Qa) and Fig.Qb). We define “weak” or “enhanced” 
pairing relative to = 0. The former is associated with 
small A or negative 1/kpa while the latter corresponds 
to either large A or large positive 1/kpa. Thus, Fig.j^c) 
is characteristic of the generic weak pairing regime while 
Fig. Ed) is characteristic of the strong pairing case pro¬ 
duced by either large 1/kpa or large A. 

We analyze the top two figures by focusing on a de¬ 
creasing abscissa which effects a transition from a trivial 
to topological phase (shown as shaded). In Fig.[^a), cor¬ 
responding to 1/kpa = 0, bz = 2.5Ep, this transition is 
driven by varying the SOC strength A. In Fig. [^b) it is 
driven directly by varying the scattering length l/kpa\ 
somewhat after the point > 0 is crossed, a further 
decrease in 1/kpa (towards the BCS limit) allows the 
system to reach a topological phase. Here we see a series 
of two transitions from topologically trivial to 4-Weyl and 
then to 2-Weyl superfluids. While there is some initial 
decline in with diminishing 1/kpa, the most signifi¬ 
cant decrease in Tf. occurs in the 2-Weyl case. 

The next two panels contrast the regime of weak pair¬ 
ing (Fig. s c)) with that of enhanced pairing (Fig.[^d)). 
In the first case, the system is BCS-like everywhere. In¬ 
creasing bz gradually suppresses and there is no clear 
signature in of the change from a trivial to a topo¬ 
logical phase (shown as shaded in the figure). As shown 
in Fig. Sd), when the pairing is enhanced, Tc becomes 
insensitive to variations in the Zeeman field until 5^ = 0. 


Shortly thereafter, the topological phase transition is 
crossed and rapidly declines. 

We can see from the last figure, in particular, that the 
satisfaction of the topological inequality and the (5/i = 
0 condition importantly define a transition (often quite 
sharp, as in Fig.j^d)) between a superfluid, characterized 
by a larger gap, and larger pair mass, Mp ~ 2m (i.e., 
more “BEC-like”), and a superfluid with a small gap, 
A/Ep <C I, and a small pair mass Mp <C m which is 
“BCS-like”. The resulting behavior of arises in the 
topological phase because there is a competition between 
the effects of a decreasing pair mass and a decreasing 
mean-field pairing gap as bz increases. The net effect is 
a lowering of in the topological phase. This can, in 
turn, be viewed as a form of BEC-BCS transition. The 
details are presented in the Supplemental Material [3T]. 

One can inquire as to why the topological transition 
becomes more apparent (as reflected in Tc) on the strong 
pairing side (Fig. [^d)), whereas it is less evident (from 
the perspective of Tc) when in the weak pairing limit 
(Fig. 0: c)). These differences are reflected in the evo¬ 
lution of the band-structure via a Van Hove singularity 
as the topological transition is crossed. To address this. 
Fig. 2 presents a constant energy contour plot for the 
band -|-T_i.k. The two axes correspond to the in-plane 
ikp) and out-of-plane {kz) momenta. For definiteness, 
we have chosen 1/kpa = 0 and /r(T), A(T) are deter¬ 
mined for a temperature just above Tc. Local extrema 
in this figure reflect Van Hove singularities, either at iso¬ 
lated points or extended in a ring-like structure. Each of 
the three panels in a given row corresponds to increasing 
values of bz with only the left-most figures in the triv¬ 
ial phase. The top three figures are in the weak pairing 
regime whereas the bottom three figures are in the regime 
of enhanced pairing. 

A key observation from these figures is that in the weak 
pairing limit there is a smooth evolution from a trivial 
to topological phase, whereas for enhanced pairing the 
band-structure evolves rather dramatically from a triv¬ 
ial and BEC-like phase to a topological and BCS-like 
phase. Indeed, the topological transition in the lower 
panel is roughly correlated with the appearance of addi¬ 
tional Van Hove singularities (as indicated). This is in 
contrast to the upper panel where Van Hove singulari¬ 
ties of the trivial and topological phases are relatively 
unchanged. These figures help interpret the behavior ob¬ 
served in Fig. 1(c) and Fig. 1(d). 

Frequency dependent spin and density response func¬ 
tions.— As in previous work m we write the correlation 
functions (above Tc) as 


(fw,q) ^ 


k ocot' 


fivEgU.) - fiv'Eg 'k+q) \ 

rjEg\^ rj F/Q'k-i-q T i^ ) 


X Waa',w'(k,k + q). 


( 6 ) 


The density-density correlation function Xpp(Q) corre- 
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FIG. 2. Evolution of the dispersion as the topological tran¬ 
sition is crossed by tuning b^. In the weak pairing limit 
(top panel), the system smoothly evolves across the transi¬ 
tion, whereas for enhanced pairing (bottom panel) there is a 
more abrupt change in band-structure. In all plots we show 
constant energy contours +E-i^]i/Ep at unitarity, with k± 
and kz in units of kp- For panels (a)-(c) we set \/kp = 0.5 
and the Zeeman field bz/Ep = 0.4,0.6,0.8, whereas for panels 
(d)-(f) we set \/kp = 1 and hzjEp = 1.2,1.7,1.8 respectively. 
Only the left-most hgures are in a trivial phase. 


spends to I = j = 0, with (Tq = 12 , whereas *, j G {x, y, z\ 
gives the corresponding spin-spin correlation function. 
The differences between the density or spin responses 
are the coherence factors (k, k-f q), which are 

rather complicated and are presented in the Supplemen¬ 
tal Material. As a numerical check on these calculations, 
the /-sum rule for the density response and related sum 
rules [H] for the spin response hold for all q. 

Quite generally, the correlation functions for a paired 
normal state can be decomposed into two parts; one in¬ 
volving the difference: q) = |A_i,k — Aii^k-i-ql 

which enters as a thermal contribution (at T 0), and 
the other involving the sum: q) = |A_i,k + 

A±i,k+q|, which we call the multiparticle contribution. 
We address the q = 0 spin response, (w, 0), (where 
i, j are x or y) so that inter-band terms dominate. Thus, 
for the ±1 subscript in the density response, the —1 band 
label yields the main contribution, whereas in the spin 
response the -1-1 band label is most important. 

Figure 3(a) shows Xs„s„(a;,0) for both the trivial and 
topological phases. In the trivial phase there are two 
clearly resolvable peaks; the first peak is associated with 
the thermal contribution and the second with the mul¬ 
tiparticle contribution. By contrast, there is only one 
peak in the topological phase. A related signature for 
the Hall conductivity (in 2D) at T = 0, rather than, as 
here, above Tc, was suggested earlier [16]. 

Importantly, this provides a means of distinguishing 
between the trivial and topological phases. We can an¬ 
alytically identify the position of the maximum in the 
first (thermal) peak, which is due to a flat band in 
^(2,-)(k, 0), and appears at precisely 2bz- The thresh¬ 
old for the second peak is wi = minkA(^’+)(k, 0). In the 


(a) Spin-spin response 



(b) Density-density response 


bJEf=2 - 



FIG. 3. Gontrast between topological (solid, red) and trivial 
(dashed, black) phases of the frequency dependent spin-spin 
(Fig. n: a)) and density-density (Fig. [^b)) correlation func¬ 
tions. Both response functions are calculated at Ijkpa = 0 
and \/kp = 1, with respective wave-vectors of q = 0 and 
q = 0.5A:_Fi for the spin and charge responses. The inset in 
Fig.j^b) shows the energy contours of Ei^’^^(k, (\)/Ep in the 
topological phase, with k±_ and kz in units oi kp. The dashed 
lines highlight the saddle point Van Hove singularity whose 
magnitude determines the frequency location of the peak re¬ 
sponse in Fig. db). 


trivial phase we find that, if /r > 0, Wi = 2A, whereas if 
^ < 0, wi = 2(A^ -I- ■ Hence wi is strictly greater 

than the frequency of the first peak (26^), thus yielding 
two distinct peaks in the response function. In the topo¬ 
logical phase, wi = 2hz so that the two peaks merge. 

We now focus on the density-density correlation func¬ 
tion Xpp(a;,q), which is only non-zero when q 0. This 
is shown in Fig. 3(b) for the case of unitarity: l/kpa = 0, 
and we can again compare the trivial and topological 
phases. Here X/kp = 1 and we plot the imaginary part 
of the response function, Xpp(q, w) [211 [31], deep in the 
topological phase {bz/Ep = 2) and in the trivial phase 
{bz/Ep = 1.2) at q = O.bkpz and T = 0.21Tp (just 
above Tc). 

In the trivial phase there are two peaks, one associ¬ 
ated with thermal contributions involving i?(^’“^(k, q) 
and the second with the multiparticle component involv¬ 
ing E^'^’E) (j^;^ q)^ jjj topological phase, there is a large 
peak at uj/Ep = 0.6, which arises from a (2D) saddle 
point Van Hove singularity contribution in A*^^’+)(k, q). 
This is associated with VkA(^’+^(k, q) = 0, which (via 
the density of states) enters as a denominator in the re¬ 
sponse functions. These saddle point Van Hove singular¬ 
ity effects are well known [23 [25 and are illustrated in 
the inset on the right and discussed in the Supplemental 
Material. Importantly, here we observe that as the sys¬ 
tem enters the topological phase they amplify the peaks 
in the density-density correlation function, thus helping 
to distinguish between the trivial and topological phases. 

Conclusions.— This paper addresses how an intrinsi¬ 
cally produced condensation temperature varies across a 
topological transition, induced by varying SOC, Zeeman 
coupling, or the scattering length. Importantly, the intro¬ 
duction of fluctuations necessarily introduces a feedback 
of the topological band-structure into Tc. The passage 
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from the trivial to the topological phase is accompanied 
by a transition in which the system is driven towards a 
low Tc, more BCS-like phase with smaller pair mass and 
smaller gap. Nevertheless, there is a range of in the 
topological phase where Tc ^ O.lTp, which is experimen¬ 
tally accessible [50] , 

We also present methods of detecting the topological 
band-structure above Tc, exploiting frequency dependent 
peaks in the density and spin responses. The topological 
transition appears in the spin response as a recombina¬ 
tion of two peaks, which are separate in the trivial phase. 
In the topological superfluid, the dynamical density re¬ 
sponse exhibits a greatly amplified peak associated with 
a (2D) saddle point Van Hove singularity. In both the 
response functions and Tc we find that the topological 
transition appears quite smooth in the weak pairing and 
much more abrupt and apparent in the strong pairing 
regime. 

Acknowledgements.— This work was supported by 
NSF-DMR-MRSEC 1420709. 

Note added.— Recently, we became aware of a comple¬ 
mentary paper that considers fluctuation effects in spin- 
orbit coupled superfluids with fixed relative population 
density using a closely related formalism [55] . 
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I. DERIVATION OF THE VERTEX FUNCTION T (Q) 

Here we present derivations of the vertex function r((3), where Q = (iuj,q), along with the coherence factor 
Vrjaa' (k, k — q) , which appear in Eq. (4) of the main text. We begin by writing the non-interacting Green’s function 
in terms of projectors as follows 


Go (K) = ^ 

Ct 


iv - ^ak’ 


(SI) 


where K = {iv, k) and ^(S(k) = (1 -I- aa^) Ul^ is a projector into the single-particle band a = ±1. The matrix C/k 

is the unitary operator that diagonalizes Hq (k) to produce the single-particle dispersion ^^k = Ck + a |h|. 

Similarly, the Nambu Green’s function Q (K) for a superfluid can be written in terms of projectors as 


G {K) = ~ '^BdG (k)] ^ 


(GJK) F{K)\_^ y^,,(k) 
\F (K) G {K)J ^ w - vEak ’ 


(S2) 


where we have used the inverse of the BdG Hamiltonian, "HEdG) to define the normal and anomalous Green’s functions 
G{K) and F{K), along with their time-reversed counterparts G{K) = icjy [G(—AT)]^ lOy and F{K) = Oy [F{—K)]^ Gy. 
The projectors Vya = are constructed from the BdG eigenvectors 



(S3) 


where Gk = diag {C/k, 14} rotates the particle (hole) sector to the spin-orbit basis with a unitary matrix C/k (14), and 
we have defined 9 = cos“^(6z/|h|), i/gk = -I- cos^ 9, and Cok = -E'Ok + Q;|h|. Note that, up to a sign, ^^k limits 

to ^Qk as A —)• 0 or —)• 0, and to Eak as A —)■ 0. 

For convenience, the 4x4 projector matrices can be expressed as four 2x2 sub-matrices as 


'P^a(k) = 


(PyaOP) (5r,a(k)\ 
yi?^a(k) S'^a(k )) ■ 


(S4) 


The Green’s function G (AT) is found from the appropriate 2x2 sub-matrix with the corresponding projector 

.. X _ ^ TT A^ok + aCk) (Aak + aryCok) aA^ sincos\ t 

’ 4AokAak V oA^ sin6»cos6» (Agk - a^k) (A^k - a?7Cak)/ ^ ^ 

We can now define a quantity X{Q), known as the pair susceptibility, which has been introduced in previous 
papers [SI]; it is a natural extension of X (0) and appears in T {Q): 


X(Q) = --Tr 


^G(A) Go {K-Q) 


K 


where Go {K) = iay [Go {—K)^iay is the time-reversed, or hole, non-interacting Green’s function. 


(S6) 
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Substituting the above definitions then gives 


X(Q) = -^Tr 


Prja (k) y^ P°, (k - q) 

iij — T) ii/ — ii, 


K r]OL 


iv - r]Ea\^ ^ iv-iu; + ^a'k-q 


= -bEE E 


2 ^ ^ iv - rjEak iv - (*W - ^a'k-q) 

k r]aoc' \ lu ^ , 


Tr[P,„ (k)P°, (k-q)] . 


Performing the summation over Matsubara frequencies reduces this expression to 


1 


^(Q) = oEE 


k rjOLOi' 


^V, + l - ivf jPak) + / (Cc'k-q)) 

P Ca'k—q) i^ 


(k, k - q). 


(S7) 


(S 8 ) 


Here the coherence factor which appears in Eq. (3) and Eq. (4) of the main text is 

v^aa' (k, k - q) = Tr [P^„ (k) P°, (k - q)] . (S9) 

The vertex function can now be defined by r((5) = [x{Q) + 5 “^] ^ . Using the expression for the susceptibility in 
Eq. (S 8 ), we obtain the vertex function r(Q) as given in Eq. (4) of the main text. The familiar gap equation can 
then also be obtained from the Thouless criterion: X(0) + g~^ = 0. 


II. THE TOPOLOGICAL PHASE DIAGRAM IN THE ENHANCED PAIRING LIMIT 

Figure SI shows the major findings concerning the topological phase diagram in the enhanced pairing limit. In 
Fig. Sl(a) the transition temperature Tc and the pairing onset or mean-field transition temperature T* are plotted as 
functions of the Zeeman field b^. Here we consider the limit of enhanced pairing through large SOC strength A, which 
leads to a more BEC-like behavior at = 0. (We could also have considered small A with sufficiently large 1/kpa 
without qualitatively affecting the result.) In Fig. Sl(b) the pair mass, which is an important component of Tc, is 
plotted as a function of b^- Finally, in Fig. Sl(c) we present a plot of the zero temperature mean-field gap as a function 
of this Zeeman field. Notable here is that for a large range of 5^, the pair mass M±{T — 0), [S2] the mean-field gap 
A(T = 0), and the transition temperature Tc, are only weakly dependent on b^. Nevertheless, Fig. Sl(a) shows a 
striking contrast between the two temperatures Tc and T*. Whereas T* decreases continuously with Zeeman field, 
Tc tends to be relatively constant until the point Sg = 0 is crossed, after which point the topological phase transition 
is encountered and Tc starts to rapidly decline. This demonstrates that T* is completely independent of topology, 
whereas Tc clearly reflects the topological phase transition. 

Now we give insight into the behavior of Tc, which depends on the solutions to the mean-field gap and number 
equations and on the behavior of the pair mass. From the perspective of Tc, we find that the trivial to topological 
transition reflects a BEC-BCS transition. This is reminiscent of other topological superfiuids, such as a Pa, -I- ipy 


(a) Phase diagram (b) Pair mass (c) Mean-field gap 



FIG. SI. Plots of the (a) Tc and T* phase diagram, (b) zero temperature pair mass, M±{T = 0), and (c) zero temperature 
mean-field gap, A(r = 0), as a function of bz in the enhanced pairing limit {l/kpa = 2 with X/kp = 1). In all figures the 
topological superfluid phase is shaded blue, and a dotted line separates the region with 5fi < 0 (left) and Sfj. > 0 (right). 
The superfluid transition temperature Tc is insensitive to bz until 5g > 0. However, the mean-field transition temperature T* 
monotonically decreases with increasing bz. The weak dependence of Tc on bz is correlated to similar behavior of the pair mass 
shown in (b) and the mean-field gap at zero temperature shown in (c). 
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superfluid, where the topologically non-trivial phase occurs on the BCS side of the BCS-BEC crossover. Ultimately 
we will show that the central physics arises from the behavior of r((5) and r(0), which, through the coherence factors 
implicitly contain the BdG wave-function. In this way the change in topology is correlated to changes in Tc- 


A. Dependence of the transition temperature Tc on the Zeeman field bz 

Our aim now is to understand how Tc depends on bz- While mean-field theory would predict a phase transition 
at T *, Tc is crucially dependent on the pair mass. To make this discussion more concrete, we refer to Fig. S2 which 
indicates how the transition temperature is deduced. The procedure is as follows. First we define the quantity Apg(r), 

^1^{T) = Y,T{Q), T<Tc, (SIO) 

Q#o 

which is commonly referred to as the pseudogap. From Fq. (5) in the main text, we see that Tc is the point at 
which this function coincides with the mean-field gap A^(T). Figure S2(a) presents a plot of the mean-field gap vs. 
temperature for two different values of which are associated with the topological phase. For these two values of b^, 
Fig. SI shows that both the pair mass and the transition temperature are decreasing rather rapidly with increasing 
bz- The dotted lines which intersect these curves are plots of Apg(r) and the points of intersection represent Tc for 
each value of bz- Also labeled is T* for each value of bz- It is also clear that the slope of the lines is greater for the 
lower bz case. This is because the pair mass is larger. Consequently the fluctuation suppression of Tc relative to T* is, 
thereby, greater. This leads to a larger separation between T* and Tc at lower bz- Nevertheless, because T* is so much 
larger here, Tc is as well. Indeed, Fig. S2(a) shows that even in the narrow range of bz, the transition temperature 
rapidly drops, as observed in Fig. SI (a). 

Figure S2(b) is a similar representation of our transition temperature calculations for two values of bz in the trivial 
phase. This regime is associated with the BFC limit where one can see that both the pair mass and Tc are very nearly 
independent of Zeeman field. This should be contrasted with the behavior of T* for these two values of bz, which 
shows a substantial decline with increasing bz- 

Of the four values of bz illustrated in Fig. S2(a) and Fig. S2(b), the higher bz = 7.6Ep is clearly in a BCS-like limit 
where the pair mass is small and Tc~T*. The lower two values of bz [bz/Ep = 2, 5) are in the BFC regime, where 
Tc and the pair mass are relatively constant. Finally, the case bz = 7.3Ep is in an intermediate regime between BCS 
and BEC-like behavior. 

We can understand these BCS and BFC endpoints semi-quantitatively. From the definition of r((5) in the text it is 
straightforward to show that aoApg(T) = (MT/27r)^/^C(3/2), where M = is the geometric mean of the 

pair mass and ({z) is the Riemann-zeta function. In the extreme BFC limit (where A(rc) « A(0), M{Tc) ~ M(0)), 
it follows that Tc ~ A^(0)/M(0), or that the transition temperature scales with the inverse pair mass. This is to be 
associated with the trivial phase. In the extreme BCS limit we have A^(Tc) = c(I — Tc/T*)A^(0) (with dimensionless 
constant c). This yields (1 — Tc/T*) = (Tc/T*)^^^ Apg(T*)/cA^(0). Solving for Tc in this small pseudogap (BCS) 
limit implies that Tc ~ T*, as expected. This describes the situation within the topological phase. Note that this 
BCS expression is independent of both the pair mass and the gap. 


B. General structure of pair mass 


The previous figures indicate that the pair mass is an important component of the present theory. The pair mass 
can be written as 


1 

2A^ 


1 9^ f — ivf (EqU.) + f (Ca'k-g)) 

2ao dqidqj ^ \ + Ca'k-q 


'^Tjaa' (k, k q) 


q=0 


(Sll) 


where Vriaa' (k, k — q) is the coherence factor defined earlier. By symmetry, the off diagonal elements {i ^ j) of 7/Mij 
vanish, and furthermore l/M^x = ^iMyy = 1/Mj_, and we also define 1/Mzz = 1/Af||. 

The most significant features of the pair mass, as shown in Fig. Sl(b), can be understood from this equation. We 
start by considering the pole structure. Of the eight denominators of the form rjEak + ^a'k-q at q —)■ 0, only the term 
with Q; = Q;' = r 7 = —1 can be zero. The critical value, be, for which this denominator can vanish is found to be /x < 0 
with &c = 5 + ItI) • Note that this condition occurs half-way between the topological phase transition 

and the Sfi = 0 transition. 
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(a) A{T) and ^g{T) in the topological phase (b) A{T) and ^pgiT) in the trivial phase 




FIG. S2. Dependence of the mean-field gap A(r) (solid curves) and pseudogap Apg(r) (dashed curves) as a function of 
temperature for the strong coupling limit {l/kpa = 2 and X/kp = 1). The left (right) panel shows the topolgical (trivial) phase 
respectively. The strength of Zeeman coupling hzjEp is indicated in the legend, and a dotted line of the same color corresponds 
to the projection of Tc where A(Tc) = Apg(Tc). Recall that by definition, T = T* defines where the mean-field gap vanishes 
(A(T*) = 0). In the left panel, the stronger Zeeman coupling, b^/Ep = 7.6, leads to a smaller T*. However, the fluctuation 
effect at h^/Ep = 7.6 is reduced and Tc is pushed toward T*. In other words, Tc « T*. In the right panel, as bz is increased 
T* drops significantly whereas Tc does not change. In this case the Apg and Tc curves almost overlap each other for both bz 
values. In fact, for this enhanced pairing limit, the pseudogap and Tc curve would be nearly identical for any bz < 6c. Note 
further that A(0) « A(Tc) in this limit. 


For < be, the system is topologically trivial. This regime can be understood if we note that no denominators in 
the inverse pair mass can be zero. Therefore, the integral is dominated by the coherence factors at small momentum, 
which will be approximately constant. Furthermore, as we will show, in this same limit the mean-field solution is also 
strongly insensitive to bz ■ This combines to produce the constant curve for the topologically trivial phase as observed 
in the pair mass, and consequently T^. 

Shortly after the topological phase transition is crossed, bz > be, a line of singularities in the integrand emerges 
defined implicitly by the curve (6^ — A^/4) = |h (k)|. After performing the q derivatives in Eq. (Sll), a term will 
appear with the pole (—£’_.k + C-.k) ^■ At larger bz, this third order singularity will dominate all other terms in the 
integrand, and the integral is well approximated by only this contribution. This will result in a rapid growth of 
at large bz, which leads to the suppression in the pair mass that was observed numerically. 

The crossover between these two regimes is characterized by a small increase in the pair mass. This can be further 
understood by noting that the lower-order poles contribute with a negative sign. This reduces the inverse pair mass 
for a small region before the third order pole begins to dominate. Thus, there is a small increase in the pair mass 
after the condition bz = be is crossed. 


C. Dependence of the mean-field equations on the topological phase transition 

The mean-field behavior of the gap function is also an important ingredient in these calculations of Te- We now 
demonstrate the sensitivity of the mean-field solutions on the condition for the topological phase transition. This 
should be contrasted with the behavior of T *, which is defined as smallest T that solves the gap equation in the A = 0 
limit: 


EE 

k oc 


(l-2/(|gka|)) 

4 ICkal 


1 + a 






1 

g' 


(S12) 


Since A does not enter into the above equation, the topological condition is not relevant for T*; this is in contrast to 
the behavior of Te- 

We now address the mean-field gap {Iq) and number (In) equations and limit our discussion here mainly to the 
ground state. This is largely for simplicity and because the topological condition is most directly reflected here. 
Moreover, it is relatively straightforward to generalize to finite T. We begin by defining functions such that the self 
consistent mean-field equations are obtained when Iq [A {bz), (bz), bz] = In [A (bz), {bz), bz] = 0. Differentiating 
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these equations with respect to b^, and applying the chain rule, one can find the change in A and /i with is 

/^\ / 3JG[A,M,bd 3JG[A,;x,fc,] \ 

I db^ ) _ _ I dA dfj, I I db^ 1 

I d/j. I I aJw[A,/j,b^] 3Jiv[A,;x,b^] I 1 OIn [A,;j,b^] I • 

\db^ / \ dA dfi / \ db^ ) 


(S13) 


We can now understand the dependence of the mean-field gap and chemical potential on the Zeeman field bz ■ Specifi¬ 
cally, at zero temperature and in the enhanced pairing limit, the gap and chemical potential are approximately indepen¬ 
dent of the strength of the Zeeman field until we cross into a topologically non-trivial phase. After this point, the gap 
begins to rapidly drop with Zeeman field, whereas the shifted chemical potential (Sp.) increases. This demonstrates 
that these changes can be associated with a BEC to BCS transition. 

Analyzing this behavior is most tractable in the limit that A —)■ 0. We consider Iq, but a similar argument applies 
for In- In this limit, the gap equation is 


/g[A {bz) ,p{bz) ,bz] 



0 



1 

~ 5 
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(S14) 


where 0{x) = 1 for x > 0 and 9{x) = 0 for x < 0 is the unit step function. These expressions show that the system 
is completely insensitive to provided that mint | \/Ck + | > bz- This condition is precisely the condition for a 

topologically trivial superfluid. As bz is increased, the system eventually crosses into a topologically non-trivial phase. 
In this phase, the step function vanishes for some range of k, allowing the gap equation to be satisfied with smaller 
A. Thus, the system becomes more BCS-like with increasing bz- 

If we consider the limit with A > 0, we see that similar arguments hold. For example, the coherence factor in Iq 
contains, in part, a term of the form 


\J (^k 4” A^) + (AfcxCk/^z) + oibz 




-I- A 2 -I- abz 


(S15) 


The a = —1 contribution can switch sign for some parameters. Specifically, the minimum condition for this term to 
change sign occurs at k± = 0. Taking this limit gives a condition that is only satisfied in a topologically non-trivial 
phase. 

In addition to the change in sign once the topological condition is crossed, if A 7 ^ 0, there is an additional variation 
with bz that is not due to the suppression of the integrand in Iq and Ijv- However, this small shift is not responsible 
for the more qualitative effect of an evolution from BEC-like to BCS-like physics as the topological phase transition 
is crossed. 


III. RESPONSE FUNCTIONS IN THE NORMAL PHASE 


In this section we derive the density-density and spin-spin correlation functions in the presence of SOC and a 
Zeeman field. Here we focus on the normal phase (T > Tc), which does not require collective-mode contributions. 
The density-density or spin-spin correlation functions, as in the main text, can be written as XsiSj (Q) = f 

(TrSqi (r) S-qj (0)), for a many-body density or spin operator Sqi = X^ss'k *^18 (o'*)ss' Ck+qs'- Here i = j = 0, with o-q = 
I 2 , corresponds to the density-density correlation function Xpp{Q), whereas i,j G {x,y,z} gives the corresponding 
spin-spin correlation function. 

We emphasize that in the normal state there is no anomalous Green’s function component, but the existence of 
normal state pairs allows one to write the correlation functions as the sum of two terms 


K 


a^G{K)aJG{K + Q) + a^F{K)(7jF{K + Q) 


Here we associate F{K) = (A/A*)_F(ir) with a pseudogap vertex contribution, which leads to 

^ ^ ^ ^ ^aa.' ,7)7)' (^7 k d" q) (ir' ^^ak) 5 

K 0.0.',7)7)' 

where we have introduced the coherence factor 


(S16) 


(S17) 


^aa',7)7)' (^5 ^ q) — (k)(7j(k “h q)] Tr(k)(7j Rt)' o' (k “h q)] • 


(S18) 
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Upon performing the summation over Matsubara frequencies, we obtain the expression for the density-density or 
spin-spin correlation function given in Eq. (6) of the main text: 


k OiO .'^ 


fivEg^u.) - /(VUa'k-rq) \ 

VEak - 1l'Ea'k+q + IW / 


'^aa' (h, k -p q). 


(S19) 


In any theory of correlation functions, it is important to ensure they satisfy the appropriate sum rules, which may 
or may not arise due to a conservation law. For the density response, the /-sum rule is a consequence of charge 
conservation, which holds even in the presence of SOC and a Zeeman field. The explicit form of the /-sum rule is 



(S20) 


Here, and below, we analytically continue the response to real frequency Xs.Sj(w, q) = lim^-^oXs-Sj (iw, q)|ii.,;=tj+i 5 , 
and then define the dissipative part XSiSj (w,q) of the correlation function. Note that for i = j, the dissipative part 
is equivalent to the imaginary part of the response function, but for i j it may be complex in general [SI]. As 
can be shown [SI], using charge conservation in the form of the Ward-Takahashi identity, the above density-density 
correlation function manifestly satisfies the /-sum rule. 

The associated sum-rule for the Xs,„Sj,(w,q) correlation function can be shown [SI] to be 

/ OO 7 

— (-ccxLsy(w,q)) = 0. (S21) 

-OO ^ 


Since XsxSy{ijJ,<l) is even in frequency, this sum rule is trivially satisfied. An additional sum rule can be derived, 
however, which is given by a different moment of the spin-spin correlation function; this sum rule is 

/ OO 7 

—X'Lsy{oj, q) = ([S'qa;, S'.qj^]) = 2i (ct^) , (S22) 

-OO ^ 

where (az) = Wz^iE)] is the spin polarization in the z direction. 

In all our numerical calculations the spin sum rule is satisfied to within one percent and similarly the /-sum rule 
for density is found to be accurate to within five percent. 


A. Van Hove singularities and indications of the topological phase 

In the main text we have noted that these two-body probes are to be contrasted with the more frequently studied 
(one-body) spectral functions, as measured in photoemission or radio frequency experiments. We consider peaks in 
the response functions as a function of the frequency uj for a fixed wave-vector q. These are associated with thermal 
contributions involving q) and multiparticle contributions involving (k, q), where we define 

A(2.-)(k,q) = |A_i,k-A±i,k+q|, 

q) = + A±i,k+q|. (S23) 

We find that the — I sign is most relevant for the density response and the -I-1 sign is most relevant for the q = 0 spin 
response. 

Van Hove singularities lead to important enhancements in the response functions because the density of states 
enters as a denominator. These singularities correspond to 

VkE(2A)(k,q) =0, (S24) 

and they are associated with minima, maxima or saddle points. Importantly, for a system with rotational symmetry, 
k± and kz form an effectively two-dimensional integrand. In such systems, it is the saddle point Van Hove singularity 
which leads to the most prominent peaks. We can now understand features of Fig. 3(b) in the main text. The 
enhanced peak of the topological phase is associated with a saddle point Van Hove singularity. The deeper one goes 
into the topological regime, the more pronounced the saddle point effects are. 

Figure S3 presents a contour plot of these constant energy surfaces for F(^’+^(k, q) which enters into the density 
response. The different panels indicate how the system evolves from the trivial to the topological phase (left to right). 
The top (a)-(c) and bottom (d)-(f) panels correspond to weak and enhanced pairing, respectively. The contours 
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FIG. S3. Evolution of Van Hove singularities (VHS) with in the contours of constant multiparticle energy (k, q). In 

the upper panel, X/kp = 0.5, and from left to right, the Zeeman field takes values hzjEp = 0.4, 0.6, and 0.8, respectively. In 
the lower panel Xjkp = 1, and from left to right bz/Ep = 1.2, 1.7, and 1.8, respectively. In both cases, the leftmost panel is 
topologically trivial, and the other two are nontrivial. Saddle point VHS are enclosed by the black dashed curves. As evidenced 
in these figures, for X/kp = 0.5, VHS can exist even before the topological phase is entered while VHS can only appear deep 
in the topological phase for X/kp = 1. In all figures, q = Q.bkpz, 1/kpa = 0, and T = 0.2iri?, and the wavevectors and 
multiparticle energy are measured in Fermi units. 


denoted by dashed lines enclose the saddle point positions. As in Fig. 2 of the main text, here we observe that in 
the enhanced pairing regime the onset of a saddle point is closely associated with the trivial to topological phase 
transition. In either case, once the system is deeper into the topological phase the saddle point Van Hove singularities 
are enhanced. This, in turn, leads to the peak enhancement in the density-density correlation function shown in 
Fig. 3(b) of the main text. 


[51] C.-T. Wu, B. M. Anderson, R. Boyack, and K. Levin, Phys. Rev. B 91, 220504(R) (2015). 

[52] The mass My is similar. 































